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Abstract 

As has been previously proposed, a minimal modification of the standard SU{3) x SU (2) x U (1) 
theory provides a viable dark matter candidate. Such a particle, a scalar gauge singlet, is natu- 
rally self-interacting — making it of particular interest given recent developments in astrophysics. 
We review this dark matter candidate, with reference to the parameter ranges currently under 
discussion. 
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I. INTRODUCTION 



Although the presence of dark matter in the universe appears to be no longer in doubt, 
the nature of the dark matter remains elusive. Spergel and Steinhardt have suggested self- 
interacting dark matter as a possible solution to numerous discrepancies between predictions 
and observations in the standard dark matter scenario. Cold dark matter predicts overly 
dense cores in galaxies, and they argue that this can be resolved if the dark matter particles 
interact with each other with a substantial scattering cross section o", while interacting only 
weakly with ordinary matter. The basic physics is easy enough to understand: collisional 
dark matter conducts heat, and as the core of a galaxy heats up it expands, thus lowering 
the central density. Furthermore, this form of dark matter may account for the discrepancy 
between theory and observation on small scale structure, the stability of galactic bars, and a 
number of other issues. In particular, Ostriker has shown that collisional dark matter 
can also account for the growth rate of massive black holes at the center of galaxies. For 
these scenarios to work it is estimated that cr/m, where m denotes the mass of the dark 
matter particle, must lie in the range 0.5-6 cm^/g, or equivalently, 2 x 10'^-3 x 10^ GeV~^. 

In this note we discuss a possible candidate within the standard SU{3) x SU{2) x U{1) 
theory. In all likelihood, the standard theory descends from a theory at higher energy 
scales, a grand unified theory, or possibly a string theory. It seems plausible that in this 
descent there would appear numerous scalar fields Xa which happen to be singlets under 
SU{?>) X SU{2) X UiX). These fields would not couple to ordinary fermions, and would be 
unknown to the standard gauge bosons. In general they can couple to Higgs fields, and 
through the Higgs fields feebly to ordinary fermions. The scalar field sector of the theory 
might read 

\ E mlxi + ^ E E r,a,xixi + ^ E E Pa,xi{v]v,l (1) 

where Lpj denotes Higgs doublets, and rjab and paj denote numerous coupling constants. The 
X^s interact with each other directly through self interaction (the 7] term), or indirectly 
through Higgs exchange (generated by the p term). Fifteen years ago, Silveira and one of 
us 0] proposed these scalar particles, Xa (called "scalar phantoms"), as candidates for the 
dark matter. We would like to re-examine this possibility in light of recent developments, 
including improvements in our knowledge of the Higgs sector. 
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II. SELF-INTERACTING DARK MATTER CANDIDATE 



Clearly the general case, with the set of parameters m^, rjab and paj, allows us a great 
deal of freedom. Following Silveira and Zee, we go to the simplest case of just one X field 
and one Higgs doublet, so that the relevant part of the Lagrangian reads 

Im'X' + \^X' + \pX\^'^^) - + \K^'^f- (2) 

We have imposed the discrete symmetry X —X, which guarantees stability of the X 
particle, as the X field is odd under this symmetry while all other fields in the universe are 
even. The coupling rj controls the self- interaction of the X, which we will choose to be in 
the range needed for cosmology. As we will see, this simplest case may be barely viable. 
After spontaneous symmetry breaking, ip acquires a vacuum expectation value v with 
= jj? /X. Writing \^{x)\ = v + H{x) as usual, we have rajj- ~ /i^ and mj^ = rn? + \p^'^- 
Replacing one of the if fields in the term pX^ by its vacuum expectation value, we 
induce a cubic coupling of the form ~ pvHX'^. Since the physical Higgs field H couples to 
quark and lepton fields, this cubic coupling leads to a coupling of X to quark and lepton 
fields. As X couples to ordinary matter only via the Higgs field, it naturally interacts feebly 
with ordinary matter and could easily have escaped detection. Indeed, at low momentum 
transfer X interacts with a quark or lepton via the exchange of the Higgs field, giving an 
interaction amplitude ~ pvf/mjj ~ pmf/m\j. As usual, / denotes generically the Yukawa 
coupling of if to a quark or lepton and rrif fv denotes the mass of the quark or lepton. 
This assumes that there is only one ip field; with multiple p fields these estimates are all 
loosened. Since, on the scale of particle physics, the masses m/ of the quarks and leptons 
in the first generation are rather small, the interaction of X with ordinary matter is further 
suppressed by a factor ~ mf/mx- In other words, as X knows about ordinary matter only 
through the Higgs, which itself couples very weakly to the first generation of fermions, its 
interaction with ordinary matter is necessarily weak, thus making X a natural candidate for 
the dark matter particle. (This also means that in the early universe, when fermions of the 
second and third generations were present in abundance, X might have played a significant 
role.) 

In its simplest version our model is governed by three unknown parameters: t], p, and 
mx- In the absence of experimental input, Silveira and Zee ||^ made the natural choice 
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that p be roughly equal to A, the Higgs self coupling, but this choice is certainly open to 
modification. 

In the model of eq. the scattering amplitude for X + X X + X receives two 
contributions: ~ rj from the direct quartic coupling, and ~ {pvY/rn^j ~ (p^yU^/A)/m^ ~ 
p^/A from Higgs exchange. The scattering amplitude is given by the larger of r] and p^/A. 
As we know nothing about 77, a natural assumption is to take r] and p^/A to be comparable, 
but again this is not required. If we take the scattering amplitude for X + X^X + Xio 
be of order ~ p^/A, then at low momentum transfer the differential cross section is given 
by a{X + X ^ X + X) ~ {p^ /Xf /m\. Inserting a/mx ~ (p^/A)^/mj^ into the Spergel- 
Steinhardt bound (2 x 10^ GeV"^ < a/mx < 3 x 10^ GeV"^) we obtain 

mx ~ 0.03-0.08 [p^/Xf'"^ GeV. (3) 

For example, for an X particle with mass comparable to the nucleon we would need Q\fX ^ 
p ^ 13-\/A. This parameter range is entirely reasonable by the standards of particle physics. 
If, on the other hand, 77 is larger than p^/A, then eq. (|^) is to be replaced by 

mx ~ 0.03-0.08 772/3 GeV. (4) 



III. COUPLING CONSTANTS 



As pointed out in the previous section, our proposed dark matter candidate in the simplest 
version is governed by three parameters: the self-coupling, 77, the mass of the particle, mx-, 
and the coupling to y?, given by p. Let us now study the limits set on their values by 
cosmological considerations. The parameters r/, p, and mx govern the temperature Ty at 
which X particles freeze out. This in turn determines the number density at freeze out, and 
thus the current number density and mass density of X particles. By demanding that a 
fraction of critical density VLx ~ 0.1 is due to X particles, we are able to estimate the mass 
of the particle. 

Following Lee & Weinberg we are to solve the evolution equation 

— = - {av}n + {av}no , (5) 

where n denotes the number density of X as a function of time t, hq the equilibrium number 
density of X, the scale size of the Universe, and {av) the thermally averaged cross section 
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for annihilation of two X particles into fermions X + X — > / + /. The first term on the 
right-hand side describes the expansion of the Universe, the second term accounts for particle 
annihilation, and the last term represents particle production. We note that the contribution 
of the cosmological constant to the critical density (estimated to be ^ 0.7 at present) 
stays constant as we go back in time, while the matter density increases as 1/ B?. At freeze 
out, therefore, the cosmological constant is negligible. 
The rate equation can be rewritten as 



df ( 45 



1/2 



mxM(/'- V), (6) 



dx {STT^fG ) 

where x = T/mx, f{x) = n/T^, fo{x) = uq/T^, G is the gravitational constant, and Nf 
is the effective number of degrees of freedom. Based on the numerical fit from Lee and 
Weinberg we obtain a numerical approximation to the value of / at freeze-out: 



f{xf) ^ 10 



45 



1/2 n -0.95 



mx{(yv) 



(7) 



The number density of X particles today is given by no = f{xf)To^, where Tq = 2.7 K = 
2.4 X 10~^^ GeV is the current temperature. The present mass density in X is given by 
Px = mxno, and so the contribution to closure density is 

Qx = f{xf)To'mx/p, (8) 

with p the present day critical density needed to close the Universe. 
The freeze out temperature, Tf = Xfirix, is given by 

1 / 45 \^/^ 

e'^''= Lr .1^ ] ^xicrv). (9) 



\lQn<^NfG 
IV. NUMBERS 

Even in the simple model of eq. (0) there is a great deal of freedom in setting values for 
the different parameters. For the sake of definiteness, we take the annihilation cross section 
{av) to be § 

lav) ^ ^-^P^ (10) 

^ ^ 47r(4mx2- 771^2)2' ^ ^ 

where mn is the mass of the Higgs, p is a coupling constant (from the X'^if'^if term in the 
Lagrangian), and is the mass of the heaviest fermion lighter than X. 



Approximating the exponent in eq. as —1, and plugging this into eq. (H), we find that 
the exphcit mass dependence drops out, B and 



ISTr^NfGTo 



Qx ~ low ^ -—— (11) 

V 45 p {av) ^ ' 



/ 

milhbarn, (12) 



h'^{(yv) 

with the Hubble constant given by Hq = 100 h km s'^Mpc^^ Utilizing our expression for 
the cross section, the contribution of X to the closure density comes out to be 

We can assume that p takes a value of order 1. For simplicity, we take p ~ ^g'^^mH/mwy, Q 
with g the weak coupling constant and mw = 80 GeV the mass of the W. For the range of 
parameters of interest to us, the heaviest appropriate fermion will be the bottom quark, so 
rrif = rrih = 4 GeV. Taking h ~ 0.7, Nf ~ 10 , and demanding that Qx ~ 0.3, gives a mass 
of 

nix ~ f^Hf^- (14) 

Plugging these values into eq. (|^), we find that the freezing temperature is given by 
Xf = Tf /m 0.04 (which depends only weakly on mx)- The particle is thus non-relativistic 
when it freezes out, as has been assumed throughout. 

It is clear that a parameter range can be found for p and mx such that the X can 
make a significant contribution to the dark matter in the Universe. However, from eqs. (^ 
and (^), our mass range implies high values for the coupling constants, suggesting that our 
perturbative approach is breaking down, and that our calculations are only to be taken 
heuristically. The simple case considered here, consisting of a single X field and a single 
Higgs doublet, is only marginally viable. As remarked earlier, the general version of our 
model allows us a great deal of freedom, and it is likely that trivial generalizations of this 
simple case could lead to viable dark matter candidates. 

Note added: In the process of writing this note we learned that similar ideas were discussed 
by M. C. Bento, O. Bertolami, R. Rosenfeld, and L. Teodoro (Phys. Rev. D 62, 041302(R) 
(2000)), and also by C. P. Burgess, M. Pospelov, and T. Veldhuis ( |hep-ph/0011335D . 
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